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Abstract 

Even if the tuning between the first and second register of a clarinet has been 
optimized by instrument makers, the lowest twelfths remain slightly too large (in- 
harmonicity) . In this article, we study the problem from two different points of view. 
First, we systematically review various physical reasons why this inharmonicity may 
take place, and the effect of different bore perturbations inserted in cylindrical in- 
struments such as bore flare, open and closed holes, taper, temperature gradients, 
visco-thermal effects, etc. Applications to a real clarinet resonator and comparisons 
with impedance measurements are then presented. A commonly accepted idea is 
that the register hole is the dominant cause for this inharmonicity: it is natural to 
expect that opening this hole will shift the position of the resonances of the in- 
strument to higher frequencies, except of course for the note for which the hole is 
exactly at the pressure node. We show that the real situation is actually more com- 
plicated because other effects, such as open holes or bore taper and bell, introduce 
resonance shifts that are comparable but with opposite sign, so that a relatively 
good overall compensation takes place. This is checked by experimental and theo- 
retical studies of the acoustical impedance curves of a clarinet. The origin of the 
observed inharmonicity in playing frequencies is therefore different, maybe related 
to the reed or the vocal tract. In a second part, we use an elementary model of 
the clarinet in order to isolate the effect of the register hole: a perfect cylindrical 
tube without closed holes. Optimization techniques are then used to calculate an 
optimum location for the register hole (without taking into account the use of the 
register hole as a B flat tone hole); the result turns out to be close to the location 
chosen by clarinet makers. Finally, attempts are made numerically to improve the 
situation by introducing small perturbations in the higher part of the cylindrical 
resonator, but no satisfactory improvement is obtained. 
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1 Introduction 



Many contributions, either experimental or theoretical, have been made in the 
past to improve our understanding of single reed woodwind instruments (see 
e.g. [1,2,3,4]). Our knowledge of the linear behaviour of the resonator is now 
satisfactory, which probably explains why most of the recent literature deals 
with the understanding of the sound production and oscillation regimes, an 
inherently non-linear problem. This does not mean that all interesting ques- 
tions concerning the linear behaviour of the resonator have been solved. For 
instance, Benade [5] proposed in the seventies some basic ideas and methods 
allowing the qualities of a wind instrument to be characterised (impedance 
peak alignment), but more detailed explorations along this line would be use- 
ful; Meynial and Kergomard [6] proposed relatively recently simple acoustical 
systems that shift the scale of a woodwind by a given micro-interval. Concern- 
ing the relation between linear properties and non-linear oscillations, ref.[7] 
discusses in general how it is possible to predict the emission frequencies, and 
even some aspects of the clean intonation and tone colour, which naturally 
leads to the design of modified instruments as soon as appropriate optimiza- 
tion criteria are defined. Recently, optimization techniques have been used in 
order to define longitudinal profiles of brass instruments [8] . 



An example of an unsolved problem is given by the clarinet. While many 
musicians agree that it is now a well tuned instrument, the twelfths corre- 
sponding to the four lowest tones remain slightly too large, by 20 or 30 cents 
[9,10,11]. One could imagine many reasons for this problem: deviation of the 
bore from purely cylindrical shape, existence of open or closed side holes (cav- 
ities), temperature gradient, etc. A generally accepted idea states nevertheless 
that the main reason resides in the register hole, because it cannot be ideally 
placed for every fingering. 



The first purpose of the present article is to analyze the origin of this tuning 
problem in detail. We will actually see that the occurrence of wide twelfths at 
either end of the clarinet register is not an obvious problem in terms of linear 
impedance theory, in contradiction with common wisdom. In real clarinets, 
various effects of opposite signs tend to compensate each other, at least par- 
tially (role of the flaring bell for instance). The second purpose of the article 
is to investigate whether or not it is possible to design a register hole allowing 
the first two complete registers to be perfectly well tuned (the definition of 
the registers will be given in the next section), ignoring all the other sources 
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for inharmonicity except the register hole. We will then reason on a simplified 
shape, cylindrical and without tone holes, the different values of the tones be- 
ing adjusted simply by choosing different lengths for the tube. Two questions 
will be discussed : 

i) what is the optimum location of the register hole for the tuning of the 
twelfths intervals between the first two registers? 

ii) is it possible to propose some simple system, located upstream from the 
highest tone hole, to provide a good correction to the residual tuning de- 
faults? 

The outline of this paper is as follows: the second section recalls some general 
features of the clarinet. The third section is devoted to the first study: it 
analyzes the length corrections and their variation with frequency for the 
main kinds of discontinuities encountered in wind instruments. It ends with 
an application to a real clarinet resonator in order to find whether the register 
hole is the main culprit of the clarinet tuning defect; calculations are compared 
with measurements. A validation of the calculations to first order in the ratio 
of the length correction to the wavelength is also given as well as a complete 
analysis of the effects of an open and closed tone holes. Then, sections 4 and 5 
are concerned with the second part of our study and answer questions i) and 
ii) respectively. Some useful formulae are given in appendix A. 



2 Generalities on the clarinet 

2. 1 Definitions 

For a given fingering, several resonance frequencies of the resonator occur, 
and several oscillation regimes can be obtained. The term register is used 
to describe the set of tones obtained for the same regime: the first register 
involves the tones corresponding to the first mode (i.e. the first resonance 
frequency) of the resonator, the frequencies being denoted fa , and the second 
register involves the tones corresponding to the second mode, the frequencies 
being denoted fa (where fa is approximately 3/i). The opening of the register 
hole allows the musician to "jump" from a given tone of the first register to 
the corresponding tone of the second register, one octave and a fifth above, the 
first two registers covering ideally more than 3 octaves 1 . For the fingerings 



We notice that while this division into two registers seems obvious to the scientist, 
the analysis of the characteristics of the tone colour is not so clear, and musicians 
divide this musical compass in four registers: chalumeau, throat, clarinet and extrem 
[9]. 
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corresponding to the higher tones of the first register, i.e. for fingerings from 
f'# to a'# 2 , the opening of the register hole is not sufficient in order to get 
the proper twelfth, and it is necessary to modify slightly the fingering in order 
to ensure the good tuning. Thus, in order to play the corresponding twelfths 
(i.e tones from c"'#), the third resonance, corresponding to an interval of a 
17th, is used, corresponding to that the physicist can call the third register 
(see figure 1). 

mode 1 I I I 



mode 2 



b' 



mode 3 



f'# a'# 



chalumeau throat clarinet extrem 



c"'# b"" 



Fig. 1. Compass and operative resonance frequency for a clarinet. 



The most commonly used clarinet is the B-flat clarinet which sounds one 
tone below the written note. Five parts can be distinguished on a clarinet: 
the mouthpiece, the barrel, the upper and lower joints and the bell. To first 
approximation, the clarinet is a cylindrical instrument, but a detailed deter- 
mination of the geometrical dimensions of the instrument reveals deviations 
from this regular, cylindrical shape (see figure 2) over parts which are short 
with respect to the wavelength. Moreover, the existence of tone and register 
holes, bore irregularities, thermal gradient and dispersion due to visco-thermal 
effects alter the resonance frequencies and their harmonicity; they must be in- 
cluded for a detailed study of the tuning of the instrument. 



2.2 Inharmonicity and sound production 



The sound of the clarinet is produced by self-sustained oscillations. Theory 
shows [13] that, at low sound intensity, the playing frequencies are determined 
by the zero of the imaginary part of the input admittance of the resonator 
(they are therefore very close to the frequencies of the maxima of the input 
impedance modulus). In addition, the reed has also a small influence on the 
playing frequency, for two different reasons: the volume velocity created by the 
motion of the reed movement, which adds to that produced by the pressure 
difference across the reed opening ; the frequency pulling effect of the reso- 
nance of the reed (at higher frequencies) combined with its damping. It can 

2 The notation system adopted in this paper is the same as that used by Baines 
[12]. 
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distance (mm) 



Fig. 2. Radius of a clarinet as a function of the distance from the reed tip (see table 
3). 

be shown that these two effects can be taken into account by introducing ap- 
propriate length corrections to the instrument, which are almost independent 
of the frequency [1,7]. At higher sound intensities, the playing frequency also 
depends on other impedance peaks; it can slightly change under the effect of 
the inharmonicity of the resonances of the resonator. A proper harmonicity of 
the first two resonances is therefore important (even for a single note); then, 
the sound frequency remains independent of the playing level, even if the first 
two resonances have similar magnitude, and good "impedance peak coopera- 
tion" [14] ensures a stable sound emission. This question is actually intricate 
because it involves the influence of many parameters. Here, for the sake of 
simplicity we limit ourselve to low intensity sounds; in other words, our goal 
will just be to achieve a perfect tuning at least for pianissimo levels. 
In this article, we will distinguish between inter-resonance inharmonicity and 
inter-register inharmonicity. The former refers to the inharmonicity between 
the first two resonance frequencies of the resonator for a given fingering; 
the latter defines the ratio between the resonance frequencies of the first 
impedance peak when the register hole is closed and the second one when 
the register hole is open 3 

3 The targeted interval is a pure twelfth, i.e. almost a ratio of 3. Actually for a 
tempered scale, the ratio is slightly different, because the tempered intervals are 
different from the natural ones: the exact value is 2 19//12 = 2.9966 . The relative 
difference is 0.11%, i.e. 2 cents (it is the difference between the tempered fifth and 
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2.3 Method and approximation 



The shape of a clarinet is very close to a pure cylinder. For a given fingering 
of the chalumeau, with first resonance frequency fi, an effective length £ e ff 
can be defined as 

4// = c/4/i , 

where c is the speed of sound. A first approximation of this length is given by 
the physical length between the input of the instrument (taking into account a 
length correction for the above mentionned reed effects) and the first open tone 
hole position, where the tube is assumed to be cut. All the perturbations from 
the cylindrical shape, such as enlargments or contractions, closed tone holes, as 
well as the precise effects of the open tone hole (depending on the downstream 
geometry and radiation property), can be regarded as small corrections to this 
first approximation leading to the correct value, £ e ff- F° r the lowest tone, for 
which all holes are closed, an equivalent length of the flaring horn is used as in 
[1]. If the resonator was purely cylindrical and lossless, the second resonance 
frequency would be exactly /2 = 3/i . However all the perturbations produce 
corrections slightly different from those on the first resonance: this creates the 
inter-resonance inharmonicity (for inter-register inharmonicity, the opening of 
the register hole also needs to be taken into account). Here, we are actually 
concerned not by the effect itself of the perturbations from the cylindrical 
shape, but by the variation of this effect between the two registers. 
The concept of frequency dependent length corrections is an adequate tool for 
this study. As in reference [6] the corrections upstream of the first open hole 
can be calculated separately for each kind of perturbation, with respect to 
the effective input of the tube (taking into account the reed effects) where the 
imaginary part of the admittance must be zero for self-sustained oscillations. 
The frequency dependence of the length corrections for all fingerings provides 
the inharmonicity. Only the effect of the part that is upstream the first open 
hole needs to be studied for each fingering. 



All calculations are carried out by ignoring the different kinds of dissipation 
(due to visco-thermal effects in the boundary layers, to radiation, etc.): weak 
dissipation is known to have a negligible effect on the resonance frequencies 
[16]. In Appendix B the effect of the resistance of a small hole is discussed, 
and even if nonlinear effects are taken into account, it is shown to be negli- 
gible. A consequence is the systematic use of purely imaginary impedances, 
which means for the resonance condition, that the input impedance is infinite. 
Perturbation to the planar mode theory is classically taken into account using 

the harmonic one, called the "skhisma"). Nevertheless this difference is very small 
and for simplicity in what follows we will ignore it (2 cents are almost inaudible, 
and within the "tunability range", especially for higher notes [15]). 
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lumped elements representing the effects of higher order, evanescent modes of 
the tubes. 



3 Length correction and inharmonicity produced by a small per- 
turbation: analysis of the different effects for a real clarinet res- 
onator 



In this section we study the analytic expressions of length corrections and 
inharmonicity of the first two resonance frequencies associated with small 
perturbations to a purely cylindrical shape. This formulation gives a good 
idea of the inharmonicity encountered on a real clarinet resonator. 



3.1 Length corrections: definition 



The effect of a geometrical perturbation on a cylindrical resonator may be ex- 
pressed conveniently in terms of a length correction to the main tube, denoted 
A£. Using an exact formulation of A£ is possible but in this study, we use an 
approximation of length corrections to first order in the ratio of the length 
correction to the wavelength: this gives a sufficiently accurate determination 
of the resonance frequencies. Moreover, the length corrections associated with 
different perturbations can be simply added as we will see in a particular ex- 
ample. The perturbation is located at a distance £ from the effective input, 
which is the origin of the coordinates once the clarinet embouchure has been 
replaced by an equivalent cylinder and the reed effects have been taken into 
account (see 3.5.1). 



3.2 Inharmonicity of the resonance frequencies 



Inharmonicity can be defined as the relative difference between the resonance 
frequency f n and n times the first resonance fi, as follows: 

TR _ fn-nf 1 _ £ eff + A£ l A£ n -A£, A£ 

nf x 4// + A4 " £ eff + [ £ eff h [ } 

where £ e ff is the acoustic length of the unperturbed system and A£ n is the 
length correction associated to the nth resonance. 

In the case IH > the basic intervals are enlarged; when IH < the intervals 
are reduced. In the present paper, results are given in cent; the cent is the 
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micro-interval equal to one hundredth of a tempered semi-tone: 



-4 



and IH) 



IH 



1 cent = 5.78 10 



cents 



5.78 10" 4 ' 



We notice that the smallest frequency deviation perceptible by the human ear, 
estimated to be 4 cents, corresponds to IH ~ 0.25%. 



3. 3 Length corrections and inharmonicity formulae for acoustic basic systems 

This section contains analysis of inharmonicity associated with the simple 
acoustical systems usually encountered in a clarinet, and depicted in table 1. 
From the calculation of first order length corrections 4 , analytical expressions 
of inharmonicity between the first and second resonance frequencies are de- 
rived. For exact formulations of length corrections, the reader is referred to 
appendix A. 

3.3.1 Insertion of a side hole 

This section uses results obtained by Meynial and Kergomard [6] concern- 
ing the effects of the insertion of an admittance branched in parallel and an 
impedance branched in series, on a cylindrical tube. Starting from the tone 
hole model proposed by Dubos et al. in [17], we propose a new calculation 
of an equivalent length that includes the effects of both the series and shunt 
impedances. Consequences for the cases of a register hole, an open side hole 
and a closed hole are subsequently discussed. 



General formulation for an admittance in parallel and an impedance 

in series We consider a side hole, or any discontinuity inserted in parallel, 
branched on a straight cylindrical tube (a typical example is the register hole 
shown in table 1). Y is the admittance of the inserted discontinuity 5 and Y up 
and Ydown are the main tube admittances upstream and downstream the dis- 
continuity respectively. The following equation can be written at the location 
of the discontinuity (x — £), 



4 The expression first order length corrections refers to length corrections expanded 
to first order in the ratio of the length correction to the wavelength. 

5 Throughout the paper the admittance is defined as a ratio of an acoustic volume 
velocity to an acoustic pressure. 
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register hole 



tone hole 



i > 



closed-hole 
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s' 


1 




— 


i 







J3 +aj2 



abrupt change in cross section 
in the upper part of the instrument 



+AJ? 



localized enlargement / contraction 
in the upper part of the instrument 



r — l 



truncated cone 



in the upper part of the instrument 




Table 1 

Acoustic basic systems and associated length corrections. 



Looking backwards to the top-end of the resonator and writing 
—jY c t&iak(£ + A£) , we get from equation (2): 



Y 



down 



-jY c tan k£ = Y- jY c tan k{i + A£) , 



(3) 
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where Ai is the length correction, Y c = S/ pc is the characteristic admittance 
of the main tube (p is the density of air and S the cross section area of the 
main tube), k = lixjjc is the wavenumber and j = . After some algebra, 
the following result is obtained: 

/ -if cos 2 H \ 
kAi = Arctan - \ Y : Y , , (4) 

V 1 - VTj m2M J 
which can be approximated in the limit of small Y j Y c by: 

kAi ~ — j — cos 2 k£ . (5) 
Y c 

For the case of a series impedance Z branched on a cylindrical tube, similar 
expressions are found to be, respectively: 



-j^sin 2 kt 
Z 



kA£ = Arctan — , (6) 

V 1 + ljf c sin2ki 

and 



fcA£~-j^sin 2 H . (7) 

Formulae (5) and (7) can also be deduced from Rayleigh's variational principle 
[18]. They are useful for different kinds of small discontinuities (see reference 
[6]), but for some particular cases, other formulae need to be derived as we 
see below. 



Side hole Referring to [17] for the tone hole model, the acoustic pressure 
and the acoustic volume velocity at both sides of the hole can be related by 
the following transfer matrix: 





where Z a and Y s are the series impedance and shunt admittance respectively. 
This formulation emphasizes the dual role of the two terms Z a and Y s , as 
opposed to formulations based on the classical T or II equivalent circuits. 
Manipulation of equation (8) gives the length correction: 

-j Y c tan k(£ + Ai) = ~ Y * + D J e , (9) 

where Ye = —j Y c tan ki . Finally: 

kAi = Arctan ( zEMJ^l^I^ ) . (10 ) 

\l + Y s Z a /A + j(Z a /Z c -Y s /Y c ) S mk£ cos k£j 1 ; 
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As a result, equation (10) is a general expression for the length correction 
of a tube branched off a cylindrical tube. It gives the influence of both the 
series and shunt impedances in the determination of the equivalent length. No 
assumption has been made concerning the dimensions and shape of the hole: 
it is therefore possible to derive first order approximations for special cases as 
tone hole, register hole or closed hole from the relative influence of Z a and Y s . 
For instance, in the case of an open hole, both Z a and Y s are inductive. With 
the low frequencies approximation, equation (10) shows that A£ is mostly 
determined by the shunt admittance. If Y s /Y c <C 1, the length correction is 
small and equation (10) gives equation (5); this is the case of the register hole. 



Open hole: register hole The register hole is a side chimney with small 
diameter but large height located a distance I ~ 140 mm down the effective 
input (typically r = 1.5 mm, h = 13 mm). Using [17] to calculate the elements 
of the model, we obtain the series impedance as: 

Z a = jup£^/S, (11) 

where v£> is the series length correction for an open hole. The shunt impedance 
is: 

Z^Y- 1 =j^(t^n[k(h + h m + h r )+jkh s }) , (12) 
Jh 

where Sh is the cross-sectional area of the hole, and h m , h s and h r are height 
corrections associated with the matching volume, the higher order modes 
impedance and radiation, respectively. Expressions for these lengths can be 
found in [17,19] (see also [20] for a correction to [17] for h s ) and are given 
hereafter: 



£(°) = -r5 2 [1.78tanh(1.84(/i + /i m )/r) + 0.940 + 0.5405 + 0.2855 2 ] _1 (13) 
h m = 1 -(l + 0.2075 3 ) (14) 

o 

h s = r(0.82 - 0.1935 - 1.095 2 + 1.275 3 - 0.715 4 ) (15) 

where 5 = r/R, k is the wavenumber of the played tone and the subscript o 
is used to refer to the open hole case. Since the low frequency limit is valid, 
the shunt impedance is well approximated by: 

Z s =j^kh>, (17) 

Oh 



where: 



h' = h + h m + h s + h r 
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This equation allows us to discuss the relative influence of the series and 
shunt impedances. Calculation gives | Y s Z a |~ 0.0045% which confirms that 
the series impedance can be ignored in the case of a long chimney with a small 
diameter. Finally, equations (10) is identical to (5), which provides: 

S h £cos 2 k£ 

kM ^~sv—> (19) 

where £ is the distance of the hole from the effective input. Equation (19) shows 
that a small open side hole always introduces a negative length correction, i.e. 
an increase in the resonance frequencies. It also shows that A£ decreases as 
the inverse of the square of frequency. Therefore, the register hole opening 
affects the first resonance frequencies much more than the second one, which 
is an important requirement for the clarinet to overblow correctly [21]. 
From equation (19), since the register hole has no effect on tones of the first 
register for which it is closed, the frequency shift of the second resonance 
frequency due to the register hole opening can be derived. We then get: 

A£ 2 2 S h cos 2 k 2 £ 
J ^-^3^/^— ' (20) 

where k 2 = 3fci, k\ being the wavenumber of the fundamental frequency of 
the played tone. This expression shows that the register hole opening gener- 
ates positive inharmonicity by pulling the second vibration mode upward in 
frequency at both ends of the register scale and has no effect at a pressure 
node as it is well known (see figure 3 where r = 1.4 mm, h = 13 mm and 
£ = 140 mm). It also appears that the two geometrical parameters which 
control inter-resonance inharmonicity are Sh'/ Sh and t. 



Open hole: tone hole In the limit of zero frequency, the shunt admittance 
of an open hole increases to infinity (see equation (5)) and the main tube be- 
haves as if it was cut at the location of the hole. The ratio Y s /Y c becomes large 
for tone holes and equation (10) can no longer be used since the perturbation 
cannot be considered as small anymore. As a consequence, A£ is no longer seen 
as an extension of the entire air column of length £ e ff but as an additional 
length to £ the distance from the effective input to the hole location. For an 
open hole of height hi (see eq.(18)), radius r at distance £ d from the open end, 
the length correction is now obtained from: 

Y 4- DY 

-jy c cotanfcAl= 'J ' , (21) 

where now Y% = —j Y c cot k£ d . Thus: 

-jY c [(l + Z a Y s /4) - jZ a Y c cot k£ d ] \ 

Y s -{l + Z a Y s /4)jY c cotk£ d ) ' 1 ' 



kA£ = Arctan 
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e f g a b c' d' e' f g' a' 
Fingering 



Fig. 3. Shift of the second resonance frequency due to a register hole of radius 
r = 1.4 mm, height h = 13 mm, located a distance i = 140 mm from the effective 
input. For the b flat fingering, the register hole does not alter the second vibration 
mode since, in this case, the register hole is located at one third the effective length. 

For a real clarinet using the expressions for Z a , Y s and Y c for an open hole, we 
see that the quantity | Z a Y s /4 | is small compared to unity, except for the four 
bottom tones for which the order of magnitude is to 0.02 . In the same way, 
the term in bracket at the numerator is close to unity. This leads us to assume 
that Z a has a negligeable effect in this study for open holes especially since we 
focus on inharmonicity. Setting the term in bracket to unity in equation (22) 
and using Taylor's formula to the third order in k£d and to the first order in 
kA£, the tone hole equivalent length becomes: 

A£ ~ 1 | , (23) 

where q = r 2 /(R 2 h') is relative to the geometry of the tone hole. From equa- 
tion (23), the length correction becomes a constant value in the lower frequency 
limit given by A£ = Idhoiel '(4 + hole) ,where £ ho i e = ti S/S h . When frequency 
increases, the tone hole length correction increases too. Then, any upper res- 
onance frequency relative to the lowest one is flattened, and negative inhar- 
monicity is generated as it can be seen in figure 4 (for a tone hole geometry 
given by r = 4.5 mm, h = 4 mm located, for each length of tube, at a distance 
£d = 30 mm from the open end). 
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Fingering 



Fig. 4. Length corrections (left) and corresponding inter-resonance inharmonicity 
(right) for a tone hole of radius r = 4.5 mm, height h = 4 mm, located at a distance 
id = 30 mm from the open end for each tube length. 



Concerning the inter-resonance inharmonicity, equation (23) leads to the 
expression 

IE ~- ~ k M .„ ■ (24) 



37r(l + g^) 2 

Discrepancy between equation (24) and numerical calculations from equation 
(22) gives less than a 4-cent difference for all fingerings except for three tones, 
d#, f and g' for which results differ of an amount to about 8 cents. This 
appears to be due to the use of the Taylor expansion of the cotan function to 
obtain equation (23). Finally, equation (24) is a simple analytic formula which 
also appears reliable, and will be used in the next. It also gives an interesting 
feature: the freedom it gives the designer to control inharmonicity due to open 
hole in the choice of the hole size. Large tone holes would be prefered to reduce 
inharmonicity of the first two resonance frequencies. 



Closed hole (or cavity) According to [17], the series and shunt impedances 
for a closed hole are 



pc 

Z s = —[-j cot k(h + h m ) + jkh s ] , 



(25) 



and 



Z a =jupt®/S, (26) 
respectively, where £ a c ^ is the series equivalent length for a closed tube given 
by 

£( c ) = - r( 5 2 [1.78cotanh(1.84(/i + h m )/r) + 0.940 + 0.540.5 + 0.2855 2 ]" 1 . (27) 
In the low frequency approximation, the shunt impedance becomes 

Z = pC ~ J 

S h k(h + hm) ' 
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and shows that it behaves as a shunt acoustic compliance C = v/pc 2 deter- 
mined by its volume v = Sh(h + h m ). Substituting Z a and Z s in equation (10) 
and denoting X = v/S£, the length correction to first order is found to be: 



A£ ~ X£ cos 2 k£ + £^ sin 2 ki ? 



(28) 



which can be rewritten in 



(X£ -£^) cos 2 k£ + £^ ] . 



(29) 



This expression shows that the effect of a closed hole inserted on a cylindrical 
tube depends mainly on its size and position. It also shows that the change in 
inertance described by £y appears in the determination of the series closed- 
hole length correction. The pratical magnitude of this term on the total closed- 
hole effect is to an amount up to 30% for the bottom holes of a clarinet; this 
cannot be neglected as it could be found in several papers [22,6]. This term 
must be taken into account especially when investigating inharmonicity of a 
real instrument for which many cavities can act for a given fingering. Since 
W<Cl, the effect of a cavity is proportional to the ratio of the inserted volume 
to the volume of air included between the reed tip and the closed hole. As a 
consequence, the virtual volume of air corresponding to the flow induced by 
the reed movement introduced by Nederveen, does not produce inharmonicity. 
With the use of equation (29), inter-resonance inharmonicity becomes: 



The result depends on two geometrical parameters, the location £ and the ratio 
v/S. Dealing with the case of real clarinet resonator and ignoring the effect of 
iffi can undersestimate the effect of cavities to about 10 cents for the lowest 
fingerings since each cavity contribution is added. This confirms the necessity 
to take into account the series closed-hole effect. Besides, as a consequence of 
the term cos 2k\£, either negative or positive inharmonicity is associated with a 
closed hole (see figure 5 for a hole volume equal to v = 0.3 cm 3 ). Finally, since 
the magnitude is proportional to the wavenumber k\, inharmonicity associated 
with a hole of fixed size and position increases with the fundamental frequency 
of the played tone. 



A£ 1 -A£ 2 = -(- 




■^ c) )(cos 2 3M-cos 2 k x £) 
£^) S m 2 2k 1 £cos2k 1 £ , 



so that: 




(30) 
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Fig. 5. Length corrections (left) and corresponding inter-resonance inharmonicity 
(right) for a closed hole whose volume v = 0.3 cm 3 as a function of k\l. 

3.3.2 Localized enlargement/ contraction 

Consider a localized enlargement (or contraction) of length £' located at dis- 
tance £ from the effective input in a cylindrical air column and let a = S'/ S 
be the ratio of the cross section area of the enlargement (or contraction) to 
the one of the main tube. Assuming a to be close to the unity (only small 
bore changes are considered), we obtain the change in resonance frequencies 
expressed through length correction according to (see appendix A eq.(A.2)): 

kA£ = (a- l)smk£'cos2k(£ + £') . (31) 

Then, the relationship between first and second resonance frequencies is: 

LH~—(a- l)(sin 3 2k x {l + £') - sin 3 2kJ) , (32) 

which states that the two diameter discontinuities located at points £ and £+£' 
create either positive or negative inharmonicity. This is shown in figure 6 in 
the case of an enlargement (a = 1.04) with a 10 cm-long inserted tube. 




Fig. 6. Length corrections (left) and corresponding inter-resonance inharmonicity 
(right) for a localized enlargement as a function of k\l (£' = 10 mm, S' = 1.04 S). 
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Substituting £ = in equations (31) and (32), the effects due to a change in 
cross section in part of a cylindrical tube are calculated. The length correction 
becomes 

kA£~ -(a- l)sin2H, (33) 



and the inharmonicity can be written 

4 



IE 



3tt 



a - 1) sin d 2k x k 



(34) 



where the term sm2k\£ still remains positive since k\i 6 [0 , 7r/2]. As a con- 
sequence, the sign of inharmonicity associated with a discontinuity depends 
only on the value of a as shown in figure 7 for a bore widening and contraction 
of 2%. 








j 









Fig. 7. Inter-resonance inharmonicity due to an abrupt change in cross section area 
as a function of k\l : contraction S' = 0.98 5 (left) and enlargement S' = 1.02 5 
(right). 



3.3.3 Change in taper close to the tube input 

The acoustical behaviour of a change in taper over a length £ c can be repre- 
sented with an equivalent electrical circuit including two inductances of op- 
posite sign and the elements of a cylindrical tube of length £ c [23]. Writing 
X\ = ljx\ and X 2 = (£ + £ c )/x2, we obtain the length correction calculated 
to first order in (i — 1, 2) (see appendix A eq.(A.3)) 



kA£ ~ X 



cos 2 k[l + 
k(£ + £ c ) 



-X, 



cos 2 kl 
k£ 



(35) 



Equation (35) states that a single taper change is equivalent to two open side 
holes with a positive and a negative inertance respectively. It is valid either 
for a positive or negative taper change, the difference being in the sign of the 
Xi which are positive for a diverging cone and negative for a converging cone. 
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Fig. 8. Geometry and symbol used for the case of a truncated cone. 

In order to evaluate the inharmonicity generated by a small truncated cone, 
it is convenient to reformulate the expression of the length correction in terms 
of two control parameters by rewriting equation (35). Since approximations 
k£ c <C 1 and £ c jx\ <C 1 are still valid, it is possible to write: 

cos 2 k(£ + 4) = cos 2 hi - k£ c sin 2k£ + o(k£ c ) , 

and 

1 1 1 

k%2 kxi{l + £ c /x\) kxi 
Therefore a simplified expression for the length correction is derived 

kA£~-— sin2H, (36) 

Xi 

where £ c /x\ and £ are the two parameters. Under these conditions and with 
the use of equation (36), inharmonicity is given by: 

IE ~ - — —sin 3 2h£ . (37) 
3n Xi 

An example of diverging cone of length £ c = 5 mm, large-end radius R = 
7.5 mm and half- angle 9 = 1.7° is shown in figure 9. Finally, looking at equa- 
tion (36), the equivalence between a positive truncated cone and an abrupt 
change in cross section for the case S' < S can be noticed. 

3.4 Other effects: radiation, dispersion and temperature 

Because of dispersion due to visco-thermal effects, the eigenfrequencies of the 
cylindrical air column cannot be exactly harmonically related. From the well 
known expression of the speed of sound with respect to frequency [24], it can 
be shown that dispersion introduces a positive inter-resonance inharmonicity 
given by 

IH = , (38) 

where r„ = ) {\f£v + (7 — l)v / 4) is the dispersion factor associated to 
the nth eigenfrequencies, £ v and £h are the viscous and thermal characteristic 
lengths and 7 the ratio of specific heats. The order of magnitude is given in 
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Fig. 9. Inter-resonance inharmonicity for a positive truncated cone as a function of 
k\l (R = 7.5 mm, 9 = IT, 4 = 5 mm). 

the next section. 

Another effect which affects the relationship between the resonance frequencies 
is the axial temperature drop. With the relation ^ = — the temperature 
gradient can be seen as a small perturbation in series which modifies the air 
density. Looking at a location x, the infinitesimal length correction is 

5x = —jZ/Z c sm 2 kx , (39) 

where Z = juApdx/S. The total length correction is obtained by an integra- 
tion over the length of the perturbed part of the tube. Assuming arbitrarily 
a linear temperature profile over the upper third of the length 4// of the in- 
strument T(x) = 3 T ° ut ~ T "> x + T in , where T in is the temperature at the input 
(x = 0) and T out is the external temperature, the effect of the thermal gradient 
is evaluated as follows: 

T i — T in r x o T in — T out r x ° 

A£ = —3 ou — / xsin 2 kx dx H — ^— — — / sin 2 kx dx , (40) 

T ou t Jo Tout Jo 

where xq = £ e ff/3. Expression (40) can be calculated analytically since the 
lower and upper bounds are given by k\£ = ir/G and kit = ir/2 and approxi- 
mation for inharmonicity is given by 

IH = -0.0422 Tout - Tin . (41) 

out 
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With these assumptions, the thermal gradient appears to produce a constant 
5-cent positive inharmonicity across the complete register. Applying a linear 
thermal gradient over the entire instrument instead of the third of its effective 
length, would lead to inharmonicity of +10 cents . Finally, with an arbitrary 
linear temperature profile which is surely not very realistic, it can be said that 
the thermal gradient have a significant effect on inharmonicity but no more 
than 5 and 10 cents over the entire scale. 

Since radiation of wind instruments depends on frequency, radiation is also a 
cause of inharmonicity but its effect is very small compared to the previous 
effects (less than 1 cent). 

3. 5 Theoretical analysis of a clarinet resonator 



Our major results concerning inharmonicity are summarized in table 2. To 



Basic perturbation 


Sign of IH 


Effective parameters 


tone hole 


< 


Sh' / Sh and Id 


closed-hole 


U < 7r/4 =>- IH > 
k£ >tt/4=> IH <0 


v/S and i 


abrupt change in cross section 


S' > S => IH > 


S'/S and £ 


in the upper part of the instrument 


S' < S => IH < 




localized enlargement / contraction 
in the upper part of the instrument 


> or < 


S'/S and £ 


diverging truncated cone 
in the upper part of the instrument 


< 


ljx\ and t 


converging truncated cone 
in the upper part of the instrument 


> 


ljx\ and I 


register hole 


> 


Sh' / Sh and I 



Table 2 



Sign of inharmonicity associated with basic acoustic perturbations close to the tube 
input. 

compare them with a real clarinet, a Buffet Crampon clarinet has been inves- 
tigated theoretically and experimentally. From the geometrical dimensions of 
the instrument, given in table 3, and with the mathematical formulations of in- 
harmonicity between the first two resonance frequencies, both inter-resonance 
and inter-register inharmonicities were predicted directly for each fingering, 
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simply adding inharmonicities associated to each perturbation. In order to 
validate these calculations, computations of the input impedance of the clar- 
inet were carried out using a transmission-line model and transfer matrices. 
Multiplying sequentially the transfer matrix of each element from the open 
end to the mouthpiece results in the input response of the instrument. From 
the obtained impedance curves, inharmonicities between the resonance fre- 
quencies were determined for all fingerings. 

Calculations with both transfer matrix and length correction approximate for- 
mulae coincide satisfactorily for the total inter-register inharmonicity as shown 
in figure 10. Focusing on the influence of a single perturbation, a good agree- 
ment is also obtained even if the two methods do not coincide for few tones. 
For a detailed validation of length correction calculations, readers are referred 
[25] . Since the two methods give similar results, our method based on adding 
corrections to first order seems appropriate. 

80 I i i i i i i i i i i i i 
60- 



40- 




e f g a be' d' e' f g' a' 

Fingering 



Fig. 10. Total inter-register inharmonicity for a Buffet Crampon clarinet. Length 
correction calculations (Q) and transfer matrix calculations (•) (the temperature 
gradient has been ignored). 

3. 5. 1 Calculation procedure and results 

The influence on inharmonicity of each bore perturbation found on a clarinet 
is now discussed and results are plotted in figures 11 and 12. Results for 
the register hole have been shown in figure 3. The principles of the length 
correction calculations are given for every kind of perturbations. 
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no 


I (mm) 


r (mm) 


h (mm) 


24 


152.7 


1.4 


13.0 


23 


166.9 


2.35 


7.1 


22 


193.4 


2.8 


7.1 


21 


202.9 


2.7 


7.1 


20 


213.2 


2.55 


7.2 


19 


230.4 


2.4 


6.4 


18 


238.9 


3.9 


10.8 


17 


241.3 


2.35 


6.8 


16 


252.2 


2.45 


8.8 


15 


270.5 


2.3 


6.9 


14 


284.3 


3.3 


8.8 


13 


287.1 


2.7 


6.8 


12 


289.3 


2.85 


6.9 


11 


308.2 


3.95 


7.4 


10 


318.8 


2.4 


7.2 


9 


348.9 


3.55 


6.5 


8 


364.8 


4.05 


8.9 


7 


369.5 


3.95 


6.1 


6 


391.0 


3.85 


8.9 


5 


413.1 


4.85 


8.9 


4 


445.6 


5.4 


6.0 


3 


472.9 


6.1 


5.0 


2 


504.9 


5.6 


5.2 


1 


543.7 


5.8 


4.4 



nAcntirm ( mm i 


U.KX111C- IC1 \l 1 VI 1 V J 


n 


n 




14 8 


117.5 


14.6 


117.5 


15.6 


152.5 


14.5 


492.5 


14.5 


526 


15.1 


565 


16.9 


578.2 


18.4 


598 


22.2 


664 


60 



Table 3 

Numerical data concerning the Buffet Crampon clarinet investigated. Tones holes 
are numbered from 1 to 24 in order of decreasing distance from mouthpiece, the 24 th 
being the register hole. Distances are given from the reed tip. For the input impedance 
measurements the mouthpiece is replaced by a cylindrical tube of length 77 mm. 
For the comparison with computations, a position correction of about —12 mm is 
introduced for the tone hole locations. 
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Fig. 11. Inter-resonance inharmonicities due to open holes (*), closed holes 
dispersion (•), temperature (+) and the flaring bell (A). 

Mouthpiece The acoustical top of the instrument (i.e the effective input) 
for measurements is defined assuming a cylinder of length Iq = 77 mm. Even if 
this can be seen as the sum of an equivalent embouchure of volume v = 11 cm 3 
plus a length correction due to reed effects equal to 10 mm as measured in 
reference [7] by comparing the resonance and playing frequencies, it must be 
noticed that this equivalent length is imposed by the experimental setup and 
ignores the variation of reed effects with frequency. This characterisation of 
the mouthpiece does not take into account the geometry of the input volume 
which influences the functioning. Actually, the theoretical analysis is made in 
two steps: 

i) the analysis of the input impedance with the equivalent embouchure (present 
section 3.5.1 ); 

ii) after a comparison of the results of the present section with experiments, 
a discussion of the effect of the conicity of the mouthpiece is given in the 
next section (see 3.6). 



Barrel and top part of the upper joint Both the barrel and the top part 
of the upper joint are tapered towards the bottom with a step discontinuity 
between them. Inharmonicity due to change in cross-section is evaluated with 
equation (34) and since the assumption X -C 1 is not true for both changes in 



23 



150 



100- A A A A a A 

A A A 

A A 

A 

| 50 

* BBSS^f! • • • * • 



.a o 

a 

o 



a -50- 



•••••••• 

* x 



**** 



-150 



• 


Barrel 


A 


Upper joint 


* 


Cross section 


-Q- 


Upper part 



e f g a be' d' e' f 
Fingering 



Fig. 12. Inter-resonance inharmonicities due to bore perturbations of the upper part: 
barrel (•), change in cross sections (*), upper joint (A) and upper part (sum of the 
three curves) (Dj. 

taper, inharmonicities are derived using exact formulae for length corrections 
(see appendix A eq.(A.3)) and equation (1). 

It appears that the barrel entails positive inharmonicity and affects much more 
the highest part of the register than the lower part. Besides, even if the change 
in cross section tends to balance the effect of the inverse conical upper joint 
(see figure 12), the total inharmonicity associated with bore perturbations at 
the input is positive and appears to be essential to correct the damage caused 
by open holes, cavities, dispersion and the flaring horn as it will be shown 
below. 

An important point to notice is that since length corrections due to taper and 
discontinuity in cross section are frequency dependent, a correction is needed 
for the location of closed tone holes, depending on the considered fingering. 



Open holes As apparent with equation (24), inharmonicity associated with 
open holes depends strongly on the length id of the tube-part below the first 
open hole. To allow realistic calculations, a good estimation of this parameter 
is nedeed. Instead of approximating id as the distance between the first two 
open tone holes, the following method is used. Starting from fingering / with 
hole 1 open, calculation of id is self-evident: an equivalent length Ai of the 
two removed admittances (Y and Yd own ) can be deduced with equation (23). 
This gives the substitution tube of length i + Ai to produce tone /. Deal- 



21 



ing with the next tone £d is now approximated by the distance between 
the first two open holes increased by the tube-piece of length A£ calculated 
previously. This allows to take into account, up to a certain extent, the effect 
of the tone hole lattice (only two or three open holes may play a role in the 
determination of the playing frequencies [26]). Then, inharmonicity is deduced 
for this fingering, a new A£ is calculated and so on. In practice, the method 
is very convenient but fails for few situations such as cross-fingering, because 
of the proximity effects between two open holes. 

As it is expected, examination of inharmonicity due to open side holes shows 
the tendency to produce negative inharmonicity across the entire scale (see 
figure 11). Besides, it appears that the highest tones of the register are more 
affected than the bottom tones because of the strong dependence of inhar- 
monicity with the wavenumber k\ (see equation (24)). 

Closed holes Inharmonicity due to a closed side-hole located above the first 
open hole is calculated with equation (30) assuming that every hole is indepen- 
dent of the others. Adding each closed hole contribution, the total closed-side 
hole inharmonicity is determined. 

Figure 11 shows the important role of the closed holes. Inharmonicity is nega- 
tive and the shape of this curve can be interpreted by studying the combined 
action of cavities located at both extremities of the tube of length £ e ff- It must 
be reminded that inharmonicity due to a closed hole is mainly proportional to 
the hole volume, the hole volumes being larger for cavities near the open-end 
than those located at the input of the instrument. On one hand, for tones 
from e to c'#, cavities located near the open-end have a larger (negative) 
effect than those located near the mouthpiece which produce small positive 
inharmonicity. As a consequence, harmonicity is more and more altered. On 
the other hand, for tones from d# to a'#, inharmonicity still remains negative 
but decreases when the chromatic scale is played which is mainly due to the 
decreasing of the number of closed holes. 

Flaring horn The expanding part at the open-end of the clarinet has been 
modeled as a catenoidal horn with a horn constant h = 0.085. Using the 
procedure given in [1] and restricting our investigation to the two lowest tones 
only, a length correction referring to the geometrical length has been assigned 
to the flaring bell according to 

A£ h = -4 + l/JfeArctan(JWi' tan(L h //i')), (42) 

where L h is the effective length including the classical end-correction 0.6Rf (Rf 
being the end-radius of the flare) and hi = h/y/1 — k 2 h 2 . Besides, a position 
correction has been given to the two holes located in the expanding part of 
the bell as Nederveen did. As expected, the flaring horn produces negative 
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inharmonicity for the lowest tone to about -15 cents and has small influence 
on fingering /. 



Dispersion and temperature Calculations of the effects of dispersion is 
obtained with equations (38) and results are plotted in figure 11. Dispersion 
effect is continuous over the entire scale of an amount to about +10 cents. 
Concerning the effect of the thermal gradient, it has been discussed in section 
3.4. 



3.6 Comparison between theory and measurements 

To correlate theoretical results and measurements, an input impedance mea- 
surement device is used (see reference [27]). It gives the acoustic linear response 
of the passive instrument (i.e without reed effects) to an harmonic excitation. 
Measurements are made in a well-insulated room and temperature is 18.5°C 
6 . The resonance frequencies are obtained by interpolating the imaginary part 
of the admittance in the vicinity of each resonance, the oscillating frequency 
satisfying 1m{Y) = 0. The absolute uncertainty on the measurements are de- 
termined to be ±5 cents. 

Figures 13 and 14 show that experiment and theory are in a good agreement, 
except for few tones for which discrepancies can be explained by an inconve- 
nient determination of the length id when investigating the open hole effect. 
Figure 14 also shows the inter-register inharmonicity with conditions close to 
the normal use of a clarinet, i.e when both the thermal gradient and the conic- 
ity of the embouchure at the end of the mouthpiece are taken into account. 
With the use of the geometry of a classical embouchure (see reference [28]), 
the effect of the conical part at the end of the embouchure has been evalu- 
ated with the exact formulae of length corrections given in appendix A (see 
eq.(A.3)). It appears that the embouchure conicity narrows the inharmonicity 
of the highest tones of the register. 

Otherwise, two clarinets were investigated experimentally. It appears that both 
clarinets have the same general tendencies in inharmonicity as figure 15 shows, 
even if small differences can be noticed for few tones. This fact has previously 
been reported in [7] and is linked to differences in bore profile of the upper part 
of the instrument. In the same reference, figure 11 shows results for two dif- 
ferent clarinets, very similar to the results of figure 15 for the inter-resonance 
inharmonicity: nevertheless a tendency for the results of the present study is to 
be 5 cents lower than the results for the clarinets measured in [7]. Moreover, 

6 Two other changes in cross section occur in the experimental setup at the in- 
put of the instrument but their effect on the resonance frequencies is substracted 
algebraically. 
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Fig. 13. Inter-resonance inharmonicity. Comparison between measurements (error 
bars) and length correction calculations (M) (uniform temperature). 
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Fig. 14. Inter-register inharmonicity. Comparison between measurements (error 
bars) and length correction calculations (M). Adding the effects of the thermal gra- 
dient and the conicity of the embouchure to (M) gives (C\). 
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adding the results for the inter-resonance inharmonicity and for the inhar- 
monicity due to register hole leads to very similar inter-register inharmonicity 
results for the two studies. 

50 1 i i i i i i i i i i i i 
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Fig. 15. Measured inter-register inharmonicities for two Buffet Crampon clarinet 
models. 



3.7 Discussion and conclusions 



i) All our calculations deals with resonance frequencies and no excitation 
mechanism has been introduced. As shown by figures 11 and 16, our analysis 
leads to the following statements concerning the origin of the general tendency 
of clarinets to have tuning deficiencies for the lowest twelfths: 

- the used linear theory does not reveal a tuning problem for the lowest tones, 

- certain dimensions of the tone holes have been chosen by makers. When 
their effects are considered as imposed, the remaining parameters to achieve 
a correct tuning are mainly in the bore profile and the location of the register 
hole, 

- a good compromise between the bell profile and the register hole location 
and size, is needed to control the tuning of the lowest twelfth (tone e), 

- the barrel and bore profile at the top end of the upper joint are of great 
importance in order to bring the lowest and middle twelfths in tune. This 
effect tends to balance the effect of the closed holes over the entire scale. 
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Fig. 16. Inharmonicities due to perturbations found in the upper part (A), due to 
hole/bell/register hole/temperature/ dispersion (o), and resulting total inter-register 
inharmonicity (sum of the two curves) (M). 

Finally, if positive inharmonicity of the resonator is observed for the lowest 
twelfths, the origin of the problem may lie in combination of the combined 
action of register hole opening, cavities, bell and inappropriate bore perturba- 
tions in the upper part. Nederveen [1] has raised the question of the necessity 
of bore perturbations in musical instruments. In the case of the clarinet, cal- 
culations reveal the need to resort to bore irregularities close to the input 
to compensate for the negative inharmonicity caused by open holes, cavities, 
flaring bell and dispersion. Reaming the upper bore profile appears to be a 
difficult task and must be done carrefully, since this part seems to control the 
accuracy of the twelfths. 



ii) As discussed in section 2, predicting playing frequencies from resonance 
curves is not an easy task. Reference [7] shows that the measured length cor- 
rection due to the reed effects (volume velocity of the reed and reed damping) 
is rather constant over the entire first register and the first half of the second 
one, then increases strongly. This allows compensation for the large increase 
of inter-register inharmonicity shown e.g. in figure 15. Nevertheless this kind 
of measurement, done using an artificial mouth, is rather difficult. It is also 
difficult to compare directly to the results obtained at different levels by an 
instrumentalist. An example of this difficulty is given in [29]. From our results 
for the two clarinets of figure 15, it is not possible to explain why the players 
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are obtaining too large twelfths for the first tones (e to g). It is a way for future 
investigation: maybe the analysis could be made using also an artificial mouth, 
but the fact is that the players often complain about too large twelfths, and 
this fact seems to be quite common: as an example, recently Buffet Crampon 
has designed a new clarinet, called Tosca, correcting the tone /. For these 
low tones, what is the optimum inharmonicity of the resonator for having well 
tuned twelfth? Is it zero? The answer is not obvious, depending also on the 
excitation level, and certainly on the "tunability range" [15,11], due to sev- 
eral factors, including the modification of both the embouchure (reed opening 
and reed damping) and the vocal tract. These musician-control parameters 
have not been extensively investigated yet. Moreover the "tunability range" 
increases with the pitch of the played tone: it is larger for the higher reg- 
ister than for the lower one, particularly due to the low level of the second 
impedance peak (for the second register, the second peak corresponds to the 
sixth peak of the first register). 

iii) In the next sections, we continue our study in the same direction but in 
a different perspective. First, the register hole is searched to allow a register 
jump for the 19-tone complete register. This is an alternative to what makers 
do since clarinettists use different fingerings to play tones from f'# to a'# and 
their associated twelfths. Second, instead of working with a real instrument, 
we consider a perfect cylindrical tube (without tone holes and flare) with only 
one register hole. The next sections are devoted to answer to the following 
general question: is it possible to find a register hole location on a cylindrical 
tube, combined with a simple perturbation in the upper part, which allows 
for the first two complete registers to be perfectly well tuned? In other words, 
we examine if it is possible to design a well-tuned instrument which is pro- 
vided with additionnal bore perturbations in the upper part which exactly 
compensate for inharmonicities due to tone holes and other effects. 



4 Optimization of the location of the register hole and tuning cor- 
rections : statement of the problem 

4-1 Assumptions 

Our goal now is to achieve harmonicity of the first two resonance frequencies 
/iand /2 , for the effective lengths corresponding to each hole, i.e. between 
two extreme values of £ e ff, written £ min and £ max - The precise intermediate 
values of the length is without importance for the present objective. As a 
consequence, the absolute effect of small discontinuities, e.g. cavities or tapers, 
is not important as well, only their relative effect between these two registers 
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being important. It is therefore convenient to consider a continuous variation 
of the length between the two extreme values. After the optimization of the 
intervals between the two registers between these two values, it will be possible 
to find the precise location of the tone holes achieving the desired scale. 



How can harmonicity between the frequency fi of the tone of the first regis- 
ter and the frequency fi of the tone when the register hole is open, be achieved 
for all lengths between £ m j n and £ max ? First we will find an optimal location 
for the register hole, located upstream of £ min , secondly we will study if a cor- 
rection system can compensate the residual defaults of this register hole. The 
dimensions of the hole are considered to be optimized by the practice of mak- 
ers, and therefore are regarded as imposed. It is actually a difficult question, 
related to nonlinear effects as well as water effects, the important fact being 
that the linear behaviour, at low level, is well known [30]. We do not take into 
account the use of the register hole as a B flat tone hole. The correction sys- 
tems are sought in order to be without manipulation by the instrumentalist, 
and therefore to act on both the first and the second registers, contrary to the 
register hole itself. 



4-2 Formulation of the optimization problem 

Optimization techniques are used to find a set of design parameters = 
can in some way be defined as optimal. These parameters 
are obtained by minimizing (or maximizing) a criterion function & which 
may be subject to constraints and/or parameter boundaries. Thus, the design 
parameters are subjected to the following requirements: 

- small changes location must be less than the distance between the reed tip 
and the first tone hole; 

- dimensions of the acoustical systems must be reasonable for the realization; 

- geometrical dimensions are positive. 

Our optimization problem is formulated as follows: 



where x is the vector of design parameters (x G M n ) and If,- and are the 
lower and upper parameter boundaries respectively, for the design parameter 




(43) 
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4-3 Criterion functions 

4-3.1 Location of the register hole 

As previously shown, the presence of the register hole means that it is prob- 
ably impossible to built a clarinet in a way that the two registers are played 
accurately in tune with the same fingerings. Thus, the objective is to find the 
register hole location that entails the smallest frequency shift of the second res- 
onance frequency for all fingerings. The first optimization is performed in order 
to play a 19-tone compass i.e a complete register. Then, a second optimiza- 
tion is performed restricting attention to fingerings from e to / ' (above this 
fingering, clarinettists do not use the same fingering to play the fundamental 
and the associated twelfth). The distance I of the hole from the effective input 
is the only design parameter. The optimization problem deals with equation 
(20) which predicts theoretically the frequency deviation associated with the 
opening of the register hole. It can be formulated with one of the two following 
criterions: 

• criterion 1: to minimize the maximum of the frequency deviation; 

• criterion 2: to minimize the mean of the square of the derivative with respect 
to k of equation (20). 

The first criterion consists in limiting the most important tuning default of 
the instrument. This criterion, which is very simple and intuitive for anyone 
who is interested in instrument design, can be written: 



An interesting point about criterion 1 is that the solution can be approximated 
analytically as it is shown in section 5.1.1. On the contrary, once the maximum 
of the deviation is achieved for a fixed register hole location, the frequency 
deviations associated with other fingerings are not taken into account for the 
evaluation of the function. 

Concerning criterion 2, the criterion function to minimize is formulated as 



where d/dk 2 is the derivative with respect to the wavenumber fc 2 , k m i n and 
k m ax being related to the playing frequency via the effective length by k min = 
3n/2l max and k max = 3ir/2£ min respectively. Contrary to criterion 1, this cri- 
terion function takes into account the deviation associated with all fingerings 
for its evaluation. The global minimum of this criterion function is achieved 
when the integrand is zero, i.e when the frequency deviation is constant for all 




(44) 




(45) 
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fingerings. Nevertheless, noting that equation (20) is necessarily equal to zero 
for a particular frequency, = 37r/2, it appears that the register hole loca- 
tion given by this criterion function is the one that minimizes inharmonicity 
variations around zero and may lead to a more homogeneous register jump 
across the entire register. 



4-3.2 Correction of the register hole effect 

The second aim of this optimization is to give suggestions on how to com- 
pensate for the register hole effect. We are looking for geometrical dimensions 
of acoustical systems whose effects alter the resonance frequencies in order to 
restore the original /2/,/1-mode frequency ratio. Similarly to the study of the 
location of the register hole and denoting IH reg and IH pert the inharmonicities 
associated with the register hole and the perturbation respectively, a criterion 
function can be written as 

criterion 1: &\ = sup (\IH reg + IH pert |) , (46) 

which deals with the maximum of the total inharmonicity. Moreover, a sec- 
ond criterion function has also been investigated by minimizing the function 
defined as 

criterion 3: = - / [IH reg + IH pert ] dk . (47) 

for which the minimum is reached when inharmonicity associated with the 
inserted perturbation is of the same magnitude as the register hole deviation 
and in the opposite direction. 



5 Results and discussions 



For the optimization, the radius of the main tube is taken as R = 7.5 mm. 
The height of the register hole is 12.5 mm and its radius 1.55 mm. The upper 
bound for the acoustical system location, is set to 154 mm which is the first 
tone hole location according to table 3. The other parameter boundaries for 
all perturbations have been chosen in order to make the realization possible. 
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5. 1 Register hole location 



5.1.1 Criterion 1: maximum of the frequency deviation 

As mentioned earlier, an optimum of criterion 1 defined by equation (44) 
can be derived analytically. In order to achieve the solution, we rewrite the 
frequency shift as follows 

and note that the magnitude of the deviation is proportional to £ and varies 
as the function F(x) = cos 2 x/x where x = k 2 l. We will first proof that the 
maximum of F(x) does not depend on £ for a certain interval of values of k 2 l 
and therefore the optimum value of IH is obtained for the minimum of £ on 
this interval. 



The effective length of the tube varies between £ max = L = 585 mm, the 
total length of the instrument including reed and radiation length corrections, 
and imin- If the compass of a register is one twelfth minus one semi-tone, the 
two extreme values of the effective length are related by £ min = £ m ax/2 18 ^ 12 
which for simplicity we assume first to be £ m i n = £ m ax/3. The location of the 
register hole, defined by £, is in the upper part of the instrument and satisfies 
£ < f ■ 

Then, for a given effective length, £ e jf lying between £ min and £ ma x, the 
wavenumber is defined as k 2 = 3n/2£ e ff, and therefore the argument of F(x) 
varies follows : 

3tt£ , „ 3tt£ 

< M < — — , (49) 



2£ 2£ 

thus 



nun 



— — < k 2 £ < — — • (50) 
2 L 2 L v ' 

Because the ratio £/L is less than 1/3, the upper bound for x = k 2 £ is therefore 

37r/2 for which F(x) = 0. Figure 17 shows the variation of F(x). From £ = L/3, 

for the interval defined by inequalities (50), F(x) varies between and 0, with 

a maximum value equal to 0.327, for x = 2.975. When £ decreases from L/3, 

figure 17 shows that the maximum values remains constant, equal to F(xq), 

except if £ becomes so small that the value of F(x) for the minimum value of 

k 2 £ reaches the same value 0.327. Using the subscripts * for referring to the 

optimal value and (i) (i = 1, 2, 3) for referring to the studied criterion, this 

corresponds to £^ = 0.205 L . Below this value, the maximum of F(x) grows 

rapidly. As a consequence, the optimum of IH for the criterion 1 is finally 

the value £+ = 0.205 L . Nevertheless, this value is not so critical, because 

the criterion function J^i for criterion 1 is linear with £, over the interval 
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Fig. 17. Perturbation function F(x) = cos 2 x/x. The same local maximum is reached 
for xq = 2.97 and x* = 0.97. 

[0.205 L, L/3] , thus the variation is not strong. Moreover, it appears clearly 
that any other criterion relative to the maximum of the frequency deviation 
will lead to a value within the interval, because the maximum is reached two 
times instead of one only. 

Finally, taking into account that the interval \l m inAmax\ corresponds actually 
to a slightly smaller interval than a twelfth i.e 2 18 / 12 , a numerical study leads 
to the optimal value for criterion 1 slightly smaller than I = 0.205 L, i.e 
= 0.2041 L . Figure 18 confirms that for criterion 1, the function J^i 
increases linearly above this optimum value but increases strongly when i 
decreases below it. 




0.05 0.1 0.15 0.2 0.25 



Fig. 18. Maximum of the frequency deviation as a function of i/L. Above 
i/L = 0.2041, the maximum value of the frequency shift increases linearly with 
the register hole location. 



The register hole deviation obtained with this "optimal" position is plotted on 
figure 19. As expected, the maximum of deviation appears at the beginning of 
the register and is almost equal to a quarter tone. Moreover, this result shows 
that the register hole is located close to its ideal position to produce the c'/g" 
transition: the register hole is located at one third of the effective length for 
this fingering so that the frequency shift is zero. 
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5.1.2 Criterion 2: mean of the square of the derivative of the deviation 



The register hole location achieved with the minimization of (45) is found to 
be 

-J- = 0.250 , (51) 

LJ 

which is not only larger than the previous result but also larger than the 
one found with numerical data given by Nederveen (£/L = 0.230). For the 
investigated clarinet (see table 3), this ratio gives £/L = 0.240. The frequency 
deviations obtained with the optimal position defined by (51) and the one 
found in literature [7,21] are very similar: the a-e" transition is correct and 
the frequency shift at the beginning of the register is about 20 cents (see figure 
19). 




e f g a be' u" e' f g' 



Fingering 



Fig. 19. Frequency deviations associated with the register hole. The register hole is 
found to be located at $p = 120.45 mm from the top end for criterion 1 (O ) and 

(2) 

l\ = 146.2 mm for criterion 2 (•). 



5.1.3 Conclusion 

Performing optimization by minimizing the criterion functions defined by (44) 
and (45) gives different results and indicates that the final location of the 
register hole is the result of a compromise. However, as expected with equation 
(48), the magnitude of the maximum of deviation increases with the distance 
I but very slightly for location larger than £/L — 0.2041. As a consequence, 
above this critical value, the location of the register hole is not so essential 
on this point of view and this can explained the difference in the register 
hole location observed between Nederveen's and the investigated clarinets for 
instance. Finally, it appears that a register hole location far from the effective 
input and above the critical value should be an interesting compromise in order 
to accurate the first twelfths of the register. Hence, we choose the location 

("2") 

given by the optimization of criterion 2 that is £\ = 146.2 mm. 
Restricting now attention to fingerings from e to /' (and not up to a'#) 
may be interesting in order to understand what instrument makers do. The 
result is that the optimization process converges to a position for the register 
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hole between the two previous extreme values for the two criterion functions. 
Moreover, for the criterion 1, the distance from the clarinet effective input is 
found to be e$P = 133.7 mm. This result corresponds very well to the location 
used by Nederveen in its calculation [1] (i.e I = 135 mm). For that case, the 
twelfths at the bottom of the scale are still very large but the maximum of the 
frequency deviation, which is also obtained for the lowest tone, has noticeably 
fallen to 35 cents (see figure 20). Concerning criterion 2, the location is very 
close from the one found for the large compass. 

50 
45 
40 

I 35 
| 30 
| 25 
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I 15 
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e f g a be' d' e' f 

Fingering 

Fig. 20. Frequency deviation associated with the register hole in the case of a real 
clarinet (i.e restricting our optimization to fingering from e to f : above f , fingerings 
are modified to ensure the well tuning and the third resonance is used). Optimal 
positions are £^ = 133.7 mm (D ) and £^ = 145.0 mm (•) 



5.2 Adjustements of natural frequencies by means of small changes of the 
bore 



5.2.1 Overview of the possibilities 

As mentioned earlier, the criterion is now to compensate for the frequency 
deviation due to the register hole by means of small changes of the bore. We 
are looking for a solution localized in the upper part of the instrument, i.e a 
solution acting for all fingerings. Looking at table 2, only three systems give 
inharmonicity in the right direction: an abrupt change in cross section area 
with S' < S, a change of conicity at the top end, and a localized enlargement 
or contraction. Concerning the case of a closed cavity, it has been shown (see 
figure 5) that both positive and negative inharmonicity can be generated: 
in order to produce negative inharmonicity, the condition k\i > ir/A which 
corresponds to playing frequencies larger than fx = c/8£ max must be valid. 
Thus, the accuracy of the twelfths at the end of the second-register scale 
(from e ' fingering) would be improved only. 
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5.2.2 Abrupt change in cross section area: S' < S 



Fixing the radius of the tube downstream the discontinuity equal to 7.5 mm, 
the upstream tube radius R! and the location I of the discontinuity are used as 
optimization variables. Performing optimization leads to the following results 



criterion 1 
criterion 3 



76.9 mm and R'^ 
76.4 mm and R'i^ 



7.1 mm. 
7.1 mm. 



The deviation from the tempered scale of the first resonance frequencies and 
inter-register for each fingering are plotted on figure 21. Results show that 
even if "optimal" parameters are quite different, the shape and the magnitude 
of final inharmonicities obtained with the two criterion functions are very 
similar. The /2//1 frequency ratio is noticeably improved for all fingerings 
except for the twelfth associated with the lowest tone which is still 20-cent 
large. An interesting result is that the discontinuity in the diameter is found 
to be located near the barrel joint which is at a distance of 60 ±4 mm from the 
closed-end on a clarinet. This result agrees with the practice of many makers: 
Nederveen noticed that change in cross section area are mostly found near the 
embouchure for reed instruments. 



-• criterion 1 




-a- criterion 3 






• • 























d' 
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Fingering 



Fig. 21. Frequency deviation for the first register (left) and inter-register inhar- 
monicity (right): abrupt change in cross section area. 



5.2.3 Localized enlargement/ contraction 

The optimization variables used in this case are the location i of the pertur- 
bation, its length £' and the radius R' as indicated on table 1. Performing op- 
timization for the case of a localized enlargement leads to the conclusion that 
small enlargement does not improve the relationship between first and sec- 
ond resonance frequencies. Even if optimization process converges with both 
criterion functions, final inharmonicities are not satisfactory. Figure 22 shows 
the tendency of the bottom notes to widen the /2//1 ratio which is musically 
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unacceptable. On the contrary, for the considered particular case, contracting 
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Fig. 22. Frequency deviation for the first register (left) and inter-register inhar- 
monicity (right): localized enlargement. 



the air column causes negative inharmonicity across the entire register and 
tends to improve the inharmonicity (see figure 23). The optimization process 
converges to the following parameters : 



criterion 1 < 



/(i) 



B! 



(i) 



62.74 mm 
-- 36.76 mm criterion 3 < 



7.17 mm 



?(3) : 

7/(3) 



E! 



(3) 



60.0 mm 
- 22.41 mm 



7.0 mm 



When minimizing the criterion function defined by equation (46), an inter- 
esting result is that the design parameters given by the optimization process 
are very similar to those of a clarinet barrel. As for an abrupt change in cross 
section, the contraction is located at a point near the barrel joint. In addi- 
tion, the acoustic length of the contraction is the same as a common barrel 
{(■barrel ~ 40 mm). When looking at what makers do, their practice varies con- 
siderably concerning the geometry of the barrel. However, some makers insert 
narrower barrel than the cylindrical portion of the entire instrument in order 
to get accurate twelfths in the upper part of the scale [31]. While the results 
given by the minimization of (47) generate negative inharmonicity (see figure 
23), the optimization process converges to two lower boundaries [R'\ = R' max 
and l\ = £ m in)- When enlarging the domain where optimization variables are 
looking for, a new minimum (lower than the one found) is achieved. Thus, we 
do not consider this solution in the next. Finally, the mean value of inhar- 
monicity is about -0.6 cents and the standard deviation of the mean value is 
equal to 11.9 cents. 
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Fig. 23. Frequency deviation for the first register (left) and inter-register inhar- 
monicity (right): localized contraction. 

5.2.4 Change in conicity at the input 

The optimization variables used for the case of tapered perturbation (see table 
1) are the radius of the upper end of our model R', the top angle 9 and 
the location point i. For each criterion function, the optimization process 
converges to the following design paremeters : 



criterion 1 < 



= 98 _ 5 mm 



3(1) 



E! 



(i) 



0.10° 

= 7.46 mm 



criterion 3 < 



ei 3) 



R! 



(3) 



73.05 mm 
0.22° 

= 7.37 mm 



The lengths of the tapered perturbations corresponding to criterion 3 and cri- 
terion 1 are about 44.5 mm and 46.5 mm which are both larger than the length 
of a classic barrel. Concerning the location point, it is near the barrel joint for 
criterion 3 and slightly below this point for criterion 1. An important thing to 
note is that the location achieved by the optimization process is always close 
to the barrel joint whatever small change is used. 

Figure 24 shows final inharmonicities obtained with the minimization of the 
two criterion functions. The /2//1 frequency ratio is improved globally across 
the entire scale for both criterion functions except for the first fingering which 
is deteriorated slightly for criterion 3. As a consequence, only the design pa- 
rameters given by the minimization of criterion 1 are considered. The mean 
value of the total inharmonicity is found to be 3.6 cents and the standard 
deviation of the mean value is about 11 cents. 

Finally, results indicate that the 20-cent deficiency still remains for the notes 
at the bottom of the scale whatever acoustic system is used in order to alter 
the resonance frequencies. Thus, no improvement can be done to compensate 
for the effect of the register hole except for tones at the end of the scale. 
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Fig. 24. Frequency deviation for the first register (left) and inter-register inhar- 
monicity (right): positive truncated cone. 

6 Conclusion 

Many deviations from the standard cylindrical tube can be observed on a 
woodwind instrument. The present work provides some theoretical results 
concerning the effects of these bore irregularities in the perspective of the 
relationship between the first two resonance frequencies for the two first reg- 
isters. The calculations using small perturbations are classical but the case of 
a side hole has been treated precisely, taking into account both the series and 
shunt impedances. 

Calculations of inharmonicity associated with small perturbations and ap- 
plication to a real case give access to the relative influence on the first two 
resonance frequencies of each pertubating system found on a woodwind and 
reveal, in the case of the clarinet, the necessity to resort to bore perturbations 
to compensate for the damage on the resonance frequencies caused mainly by 
the presence of closed holes and the register hole. 

In addition, the use of simple optimization techniques has been shown to be 
useful in order to give a more general view of the problem. However, opti- 
mization is very sensitive to the definitions of the criteria; for our purpose, 
performing optimization on the accuracy of the twelfths does not impose a 
unique possible criterion. For instance, the use of optimization techniques 
with weight for the lowest tones is of course possible and would probably give 
other results. 

The main results of the present work can be summarized as follows: 

- the origin of the occurence of wide twelfths for the lowest tones of a clarinet 
is not to be found in shifts of the resonance frequencies of the resonator, 
as opposed to common wisdom. A way for future investigations may be the 
influences of the blowing parameters, reed parameters and the vocal tract 
since players can modify both of them and claim to do it. 

- with a single corrective system, designing a new instrument having for the 
two registers a compass of a twelfth with the same fingering except the 
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register hole opening, is impossible without tuning deficiencies; 
- if the aim is the improvement of a common clarinet, the use of contractions 
near the barrel can help to some extent; nevertheless, no simple solution 
located in the upper part of the instrument can ensure the perfect tuning 
of the twelfths corresponding to the first lowest tones. 

Finally, combining global and local solutions, i.e solutions acting mainly on 
one particular tone, is necessary in order to achieve the best tuning. 
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Appendix 

A Exact calculation of length corrections 

The aim of the appendix A is to give some precise definitions and exact formu- 
lae for some elements often encountered on wind instruments. Exact formulae 
can be useful for the computation of resonance frequencies by using an iter- 
ative procedure. The geometries of the different elements are shown in table 
1. 

Abrupt change in cross section 

Looking for a equivalent system of section S, one can write at the discontinuity 

-jY' c tan ki = -jY c tan k(£ + A£) 

where Y' c = S' / pc and Y c = S/pc. Manipulating the above expression gives 
result in 



where a = S' / S is the ratio of cross sections of the tubes located upstream 
and downstream the discontinuity. 




(A.l) 
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Localized enlargement/ contraction at the input of the instrument 

Considering a localized enlargement or contraction of length £' and cross sec- 
tion area 5" located at point £ and denoting a = S' / S the change in diameter 
where S is the cross section area of the main tube, it can be written: 

Y upi = —jY c t&nk£ at point £ 
Ydown 2 = -jY c tan k(£ + £' + A£) at point £ + £' . 

which results in after some algebra 
kA£ = Arctan (tan k£' 



-T tan k£' + a (1 + TT 2 ) + a 2 tan k£'T 2/ 



where T = tan k£ and T 2 = tan k(£ + £'). 

Denoting a = 1 + X, the expression of the length correction becomes 

7A „ . ( v . .„ cos k(2£ + £')+X cos hi cos k(£ + £') 

kA£ = Arctan X sm k£ 



1 + X cos 2 k£' + X(X + 2) sin k£' cos k£ sink(£ + £') ' 

(A.2) 

In the limit that a is close to the unity, an approximation of equation (A.2) 
is given by equation (31). 



Truncated cone at the input 

This paragraph deals with the case of a perturbation in taper of length £ c 
located at distance £ from the reed tip. Calculations for a positive and a 
negative taper are identical. The unique difference is that the sign of both 
quantities X and X' has to be inverted. 

Assuming the quantities to be positive and x 2 > x±, and noting that the two 
inductances located at points £ and £+£ c are positive and negative respectively, 
the following system is derived : 

Y up2 = Y 2 + Y down2 at point £ + £ c (x 2 > 0) and L > , 
< Y upi = —Yi + Y down at point £ (x\ > 0) and L < , 
Y downl = -jY c tan k(£ + £ c + A£) , 

where A£ is the global length correction. 

£ £ + £ c 
Denoting X = — and X' = -, the result is 

Xi x 2 

kA£ = Arctan A , (A. 3) 
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where 



A 



-y^ cos 2 k(£+£ c ) yi cos 2 fcl , y yl cos ki cos k(£+£ c ) • -yn 

_ k(e+e c ) _ kt -r^^ M k(e+e c ) ^^^c 

\X'^ + iX^gtM + XX'^^A tan k{£ + Q sin U c 



In the limit of small perturbation, the length correction for a truncated cone 
is the sum of the length correction associated with both changes in taper. 



B Effect of the resistive term in the hole impedance on resonance 
frequencies 



The aim of this appendix is to give an analytic expression of the playing 
frequency in term of the resistive part R of the hole impedance. This can be 
interesting especially when resistive effects are large and can occur at high level 
when nonlinear effects appear, especially for narrow holes, like the register 
hole. This appendix describes the calculations which lead to the following 
expression for the frequency of the played tone 

/= — A ! c (B.l) 
J 1+A v ; 

where f a and f c are the frequencies obtained when the resistance term is zero 
and infinite respectively (when R is infinite, the effect of the hole disappears 
as it was closed), and A is a dimensionless factor given by 



.4 



/2R\ 2 1 Li 



V n J i ho i e e d + A£ 



with L the total length of the resonator, I the hole location, Ihoie = h'S/Sh 
and A£ the equivalent length of the hole. As done in this study for the case 
of open hole, the series impedances Z a are assumed to be small (see section 
3.3.1). 



Using dimensionless impedance quantities, the impedance equivalent Z eq to 
the tone hole impedance and to the tube downstream the discontinuity can 
be written at the discontinuity point as 

where Z h = R + jX is the tone hole impedance and Z down the impedance of 
the tube donwstream the discontinuity. Then, the input impedance Zj„ of the 
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system can be derived as follows: 



2 • 



-t + jZ h T 



l+t 

T-t 



(B.2) 



T-t 



+ J 



where 

• Zi n is the input impedance of the system ; 

• Zh = R + jX is the tone hole impedance ; 

• t — tan k£ ; 

• T = tan kL . 

For the resonance frequencies, the imaginary part of the input impedance 
vanishes, i.e 

Tm(Z in ) = + (*^^ + XT) +TR 2 = 0. (B.3) 

Denoting 

A = sin k(L - 2£) + sin kL + 2x cos kL , 
.B = cos k(L — 2^) — cos kL + 2% sin kL , 
C = 2i? 2 sin2A;L , 
where x — kihoie, equation (B.3) can be rewritten in the form 

AB + C = . (B.4) 

By writing kL = k\L + ei where k\ is the wavenumber of the fundamental 
frequency of the played tone defined as A = (when R vanishes) : 

sin kx(L - 21) + sin k x L + 2\ cos k x L = , (B.5) 

the use of Taylor's formula to the first order for the A and B quantities results 



m: 



A = €i[(L — 21) cos k\[L — 2i)+L cos k x L + 2i hoU cos k x L - 2\L sin k x L\ , 

(B.6) 



B = cosk x (L — 21) — coskiL + 2xsin k x L 

+ ei [L sin k x L - (L - 21) sin £4 (L - 2£) + 24 oZe sin k x L 

+ 2 X L cos k x L]. (B.7) 
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Neglecting terms of second order in e\ and using equation (B.5), the quantity 
AB becomes: 



AB = ( — — [L(l + cos 2k x t) + 2£ hole cos 2 kL - 21 cos k(L - 21) cos kL 

\COS klihole j 



cos 2k\k 



cos k\L 



(B.8) 



id A£ 

In the limit of k\£& <C 1 (i.e — and are smaller than unity) and noting 

kl = WTAfy we have: 



cos " lL = "2£(4^T£) +o( 7 ) 

cos2M«i=l + o(j) 

cosfci(L-2£) = sin- = + o(y) 

cos2M= -1 + • 
Substituting these results in equation (B.4) leads to the following equation 

\ k 2 thole d ~— n — ) ei + 2i? 2 sin 2kL = . 



IT 

Noting fcL = + e2 where fc c L = — (& c is the wavenumber when the 
hole is closed), it can be written : 

^hoie^—ik - h) + AR 2 L(k -k c ) = 0, (B.9) 



therefore, 



where 



/ = !j y^ , (B.io) 

MvJw^at (R11) 

/ is the frequency of the tone with no resistance, and f c is the frequency of 
the tone with an infinite resistance term. 

A careful limiting process shows that all is in order in the two extreme cases 
i.e when R — > and R — > oo. In the limit that R goes to zero, equation (B.IO) 
gives the frequency to tend to f a , the frequency of the tone when the hole is 
open. In the opposite limit when R tends to infinity, equation (B.l) gives the 
frequency to be f c . Figure B.l shows that the playing frequencies given by 
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equation (B.10) coincide quite well with the playing frequencies obtained by 
the zero values of the imaginary part of the input admittance. Moreover, in 
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Fig. B.l. Playing frequency as a function of the real part of the tone hole impedance 
(R = 10 mm, r = 7 mm,h = 9 mm). In the limit that R — > 0, the frequency of the 
played tone tends to f Q , the frequency obtained when the hole is open. In the limit 
that R — > oo, the frequency of the played tone approaches f c , the frequency when the 
hole being closed. 

the limit that R is much smaller than unity, it follows 

f = f + A(f f -f ). (B.12) 

With the use of experimental values obtained by Dalmont et al. [20] about 
the non linear behavior of an open side hole, it appears that the resistive term 
of the shunt impedance has small effects in the determination of the playing 
frequency. For instance, a 5-cent difference is obtained with R/Z c = and 
R/Z c = 0.05. 
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